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Abstract 

By using some techniques such as spectral distribution and stratification associated 
with the graphs, employed in [1, 2] for the purpose of Perfect state transfer (PST) of a 
single qubit over antipodes of distance-regular spin networks and PST of a d-level quan- 
tum state over antipodes of pseudo-distance regular networks, PST of an m-qubit GHZ 
state is investigated. To do so, we employ the particular distance-regular networks (called 
Johnson networks) J (2m, m) to transfer an m-qubit GHZ state initially prepared in an 
arbitrary node of the network (called the reference node) to the corresponding antipode, 
perfectly. 

Keywords: Perfect state transferenc, GHZ states, Johnson network, Strat- 
ification, Spectral distribution 
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In the interior of quantum computers good communication between different parts of tlie 
system is essential. The need is thus to transfer quantum states and generate entanglement 
between different regions contained within the system. In quantum information processing 
(QIP) protocols, quantum spin systems [3, 4, 5] (or quasi-spin [6]) systems serve as quantum 
channels. The idea to use quantum spin chains for short distance quantum communication 
was put forward by Bose [7] . He showed that an array of spins (or spin like two level systems) 
with isotropic Heisenberg interaction is suitable for quantum state transfer. In particular, 
spin chains can be used as transmission lines for quantum states without the need to have 
controllable coupling constants between the qubits or complicated gating schemes to achieve 
high transfer fidehty. After the work of Bose, the use of spin chains [8]- [23] and harmonic chains 
[24] as quantum wires have been proposed. In the previous work [1], the so called distance- 
regular graphs have been considered as spin networks (in the sense that with each vertex of a 
distance-regular graph a qubit or a spin was associated ) and perfect state transfer (PST) of a 
single qubit state over antipodes of these networks has been investigated, whereas in Ref.[2] the 
authors have been considered the PST of a d-level quantum state over antipodes of pseudo- 
distance regular graphs. In both of these works, a procedure for finding suitable coupling 
constants in some particular spin Hamiltonians has been given so that perfect transfer of a 
quantum state between antipodes of the networks can be achieved. The present work focuses on 
the PST of the m qubit maximally entangled GHZ states \'iIj)ghz = T^d 00_^:^) + I il;;^))- 

m rn 

To this end, we will consider the Johnson networks J(2m, m) (which are distance-regular) as 
spin networks. Then, we use the algebraic properties of these networks in order to find suitable 
coupling constants in some particular spin Hamiltonians so that perfect transfer of m-qubit 
GHZ states between antipodes of the networks can be achieved. 

The organization of the paper is as follows: In section 2, we review some preliminary facts 
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about graphs, their stratifications , spectral distribution associated with them and underlying 
networks derived from symmetric group Sn called also Johnson networks. Section 3 is devoted 
to perfect transfer of m-qubit GHZ states over antipodes of the networks J(2m,m), where a 
method for finding suitable coupling constants in particular spin Hamiltonians so that PST be 
possible, is given. The paper is ended with a brief conclusion. 

2 Preliminaries 

In this section we recall some preliminary facts about graphs, their stratifications, spectral 
distribution associated with them and underlying networks derived from symmetric group Sn 
called also Johnson networks. 

2.1 Graphs and their stratifications 

A graph is a pair F — {V,E), where F is a non-empty set called the vertex set and is a 
subset of {{x,y) : x,y & V,x ^ y} called the edge set of the graph. Two vertices x,y & V are 
called adjacent if (x, y) e E, and in that case we write x ~ y. For a graph F = {V,E), the 
adjacency matrix A is defined as 



Conversely, for a non-empty set V, a graph structure is uniquely determined by such a matrix 
indexed by V. 

The degree or valency of a vertex x &V is defined by 



where, | • | denotes the cardinality. The graph is called regular if the degree of all of the vertices 
be the same. In this paper, we will assume that graphs under discussion are regular. A finite 
sequence Xo,Xi, ...,Xn G V is called a walk of length n (or of n steps) if Xi^i ~ Xj for all 






otherwise 



i^i^) ^ \{y &V : y x}\ 



(2-2) 
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i = 1,2, ...,n. Let 1^{V) denote the Hilbert space of C-valued square-summable functions on 
V. With each P e V we associate a vector such that the /3-th entry of it is 1 and all of 
the other entries of it are zero. Then {\/3) : (3 E V} becomes a complete orthonormal basis of 
l'^{V). The adjacency matrix is considered as an operator acting in 1^{V) in such a way that 

m = i: (2-3) 

Now, we recall the notion of stratification for a given graph F. To this end, let d{x, y) 
be the length of the shortest walk connecting x and y for x ^ y. By definition d{x, x) = 
for all X E V. The graph becomes a metric space with the distance function d. Note that 
d{x, y) — 1 \l and only if x ~ We fix a vertex o e F as an origin of the graph, called the 
reference vertex. Then, the graph V is stratified into a disjoint union of strata (with respect 
to the reference vertex o) as 

oo 

V=[j r,(o), r,(o) ■.= {aEV: d{a, o) = t} (2-4) 

i=0 

Note that Ti{o) = may occur for some i > 1. In that case we have Ti{o) = Ti^i{o) — ... — 0. 
With each stratum Ti{o) we associate a unit vector in l'^{V) defined by 

where, = |ri(o)| is called the i-th valency of the graph (ni :— [{7 : 9(o, 7) = i}\ — |rj(o)|). 

One should notice that, for distance regular graphs, the above stratification is independent 
of the choice of reference vertex and the vectors \(f)i),i — 0,l,...,d — 1 form an orthonormal 
basis for the so called Krylov subspace Kd{\(l)o),A) defined as 

Xd(|0o),A) =span{|0o),A|0o),---,A^-Vo)}. (2-6) 
Then it can be shown that [25], the orthonormal basis \4>i) are written as 

|0,) = P,(A)|0o), (2-7) 

where Pj = ao + aiA + ... + aiA^ is a polynomial of degree i in indeterminate A (for more details 
see for example [25, 26]). 
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2.2 Spectral distribution associated with the graphs 

Now, we recall some preliminary facts about spectral techniques used in the paper, where more 
details have been given in Refs. [26, 27, 28, 29]. 

Actually the spectral analysis of operators is an important issue in quantum mechan- 
ics, operator theory and mathematical physics [30, 31]. As an example ii{dx) = \ip{x)\'^dx 
{lJ>{dp) — |'0(p)po?p) is a spectral distribution which is assigned to the position (momentum) 
operator X{P). Moreover, in general quasi-distributions are the assigned spectral distribu- 
tions of two hermitian non-commuting operators with a prescribed ordering. For example the 
Wigner distribution in phase space is the assigned spectral distribution for two non-commuting 
operators X (shift operator) and P (momentum operator) with Wyle-ordering among them 
[32, 33]. It is well known that, for any pair (^4, \4>o)) of a matrix A and a vector \4>o), one can 
assign a measure /i as follows 

fi{x) = {cl)o\E{x)\(Po), (2-8) 
where E{x) — J2i \ui){ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = J xE{x)dx. (2-9) 

Then, for any polynomial P{A) we have 

P{A) = j P{x)E(x)dx, (2-10) 

where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-8) and (2-10), the expectation value of powers of adjacency matrix A over 
reference vector |0o) can be written as 

(^ol^^'l^o) = / xXdx), m = 0,l,2,.... (2-11) 

Obviously, the relation (2-11) implies an isomorphism from the Hilbert space of the stratifi- 
cation onto the closed linear span of the orthogonal polynomials with respect to the measure 
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From orthonormality of the unit vectors \4>i) given in Eq.(2-5) (with \4>q) as unit vector 
assigned to the reference node) we have 

5ij = (0,10,) = / Pi{x)Pj{x)ii(dx). (2-12) 

By rescahng Pk bs Qk — y/^^i ■ ■ -^^kPk, the spectral distribution /i under question will be 
characterized by the property of orthonormal polynomials {Qk} defined recurrently by 

Qo{x) = 1, Qi{x) = X, 

xQkix) = Qk+i{x) + akQk{x) + u}kQk-i{x), k>l. (2-13) 
The parameters and cuk appearing in (2-13) are defined by 

ao = 0, K-bk- c^, = Pi ^ h-iCk, /c = 1, d, (2-14) 

where, k = k\\s the degree of the networks and fej's and Cj's are the corresponding intersection 
numbers. Following Ref. [34], we will refer to the parameters and uj^ as QD (Quantum 
Decomposition) parameters (see Refs. [26, 27, 28, 34] for more details). If such a spectral 
distribution is unique, the spectral distribution ^ is determined by the identity 

^''^^^ Irx-v X - an - ^ ^1 ^2 ^ Qd+i{x) ^x-x^' 

X — OLO ^ 

^ X — Oi'^ 

where, xi are the roots of the polynomial Qd+i{x). Gfj,{x) is called the Stieltjes/Hilbert trans- 
form of spectral distribution ji and polynomials {Q^^^} are defined recurrently as 

g?)(x) = l, Q^^\x)^x-a^, 

xQ^^\x) = Q^l,{x)+ak+iQ^\x)+uj,+^Q^Ux), k > 1, (2-16) 
respectively. The coefficients 7; appearing in (2-15) are calculated as 

7,:= \im[(x-xi)Gi,{x)] (2-17) 



PST of m-qubit GHZ states 



8 



Now let G^{z) is known, then the spectral distribution fi can be determined in terms of 
X;, Z = 1, 2, ... and Gauss quadrature constants 7;, Z = 1, 2, ... as 



(2-18) 



1=0 



(for more details see Refs. [35, 36, 37, 38]). 



2.3 Underlying networks derived from symmetric group Sr, 



Let A = (Ai,...,Am) be a partition of n, i.e., Ai + ... + \m = n. Wc consider the subgroup 

Sm ® Sn-m of Sn with m < [^1. Then we assume the finite set = -5—^^ with \M^\ = 

m\(n-m)\ ^ vertex set. In fact, is the set of (m — l)-faces of (n — l)-simplex (recall that. 



the graph of an (n — l)-simplex is the complete graph with n vertices denoted by Kn). If we 
denote the vertex i by m-tuple (^1,^2, ■■■,^m)) then the adjacency matrices A^, A; = 0, 1, ...,m 
are defined as 



(A 



, A; = 0, 1, m. 



(2-19) 



1 if 9(i,j)=k, 
otherwise (i,j e M-^) 

where, we mean by d{i,j) the number of components that i — (ii, ...,im) and j = (ji, ...,jm) 
are different (this is the same as Hamming distance which is defined in coding theory). The 
network with adjacency matrices defined by (2-19) is known also as the Johnson network 
J{n, m) and has m-\-l strata such that 



/ 



Kq = 1, Ki 



m 
m — I 



n — m 



V 



, /=l,2,...,m. 



(2-20) 



One should notice that for the purpose of PST, we must have Km — ^ which is fulfilled if 
n — 2m, so we will consider the network J(2m, m) in order to transfer m-qubit GHZ states 
(hereafter we will take n = 2m so that we have = !)• If we stratify the network J(2m, m) 
with respect to a given reference node |0o) = \h,i2, ■■■,im), the unit vectors i — 1, ...,m 
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are defined as 



m 

\h) ^ J2 \ii,--ii-i,i'i,ii+i,-;ik-i,i'k,^k+i--;ir 



I0„) = -^ E K,4-,4)- (2-21) 



Since the network J(n, m) is distance-regular, the above stratification is independent of the 
choice of reference node. The intersection array of the network is given by 

bi = {m-lf ; ci^P. (2-22) 

Then, by using the Eq. (2-14), the QD parameters ai and Ui are obtained as follows 

ai^2l{m-l) 1^0,1,. ..,m; u;i^f{m-l + lf ,1^1,2,. ..,m. (2-23) 

Then, one can show that [27] 

A\(f)i) = (/ + l)(m - 010m) + 2/(m - l)\cf>i) + l{m-l + l)|</./_i). (2-24) 

3 Perfect Transfer of m-qubit GHZ states over antipodes 
of the network J(2m, m) 

The quantum state transfer protocol involves two steps: initialization and evolution. First, a 
m-qubit GHZ state 1'^)^ = :^(|00- ■ -0)^ + |11 . ■ . 1)a) G 'Ha to be transmitted is created. 



The state of the entire system after this step is given by 

1 

71 



1^(^ = 0)) = |V)a|00.^)b = 4^(|00...0)^|00...0)b + |11...1)a|00...0b)). (3-25) 
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Then, the network couphngs are switched on and the whole system is allowed to evolve under 
U{t) — e""'^* for a fixed time interval, say to. The final state becomes 

mo)) = ^(|00 . . . 0)^|00 . . . 0)b + E fjAih) \j)) (3-26) 

where, \j) = |jij2---Jm) = |0 . . . 0^0 . . . 0^0 . . . 0^0) and \A) = \ n^00_^) so 

that fjA{to) '■= (j|e~'^*°|A). Any site S is in a mixed state if \ fAB{to)\ < 1, which also implies 
that the state transfer from site ^4 to 5 is imperfect. In this paper, we will focus only on PST. 
This means that we consider the condition 

|/AB(io)| = l for some < to < oo (3-27) 

which can be interpreted as the signature of perfect communication (or PST) between A and 
B in time t^ (with \B) = |00 ... Oil ... 1)). The effect of the modulus in (3-27) is that the 
state at 5, after transmission, will no longer be |^), but will be of the form 

^(|00...0) + e*^|ll...l)). (3-28) 
v2 

The phase factor e"^ is not a problem because (f) can be corrected for with an appropriate phase 
gate (for more details see for example [9, 16, 39, 40]). 

As regards the arguments of subsection 2.2, the evolution with the adjacency matrix H — 
A = Ai for distance- regular networks starting in |0o); always remains in the stratification 
space. For distance-regular network J {2171,171) for which the last stratum, i.e., \(f)m) contains 
only one site, then PST between the antipodes |0o) and \(j)m) is allowed. Thus, we can restrict 
our attention to the stratification space for the purpose of PST from |0o) to 10^)- 

The model we will consider is the network J {2m, m) consisting oi N — C?I^ — sites 

\ ' / o m m!m! 

labeled by {1, 2, A^} and diameter m. Then we stratify the network with respect to a chosen 
reference site, say 1. At time t = 0, the state in the first (input) site of the network is prepared 
in the m-qubit GHZ state |V'm)- We wish to transfer the state to the A^th (output) site of the 
network with unit efficiency after a well-defined period of time. We shall assume that initially 
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the network is in the state |00 . . . 0)^|00 . . . 0)^. Then, we consider the dynamics of the system 
to be governed by the quantum-mechanical Hamiltonian 



m 



k=l l<i<j<2m 



(3-29) 



where, (7j is a vector with famihar Pauh matrices af, erf and af as its components acting on 
the one-site Hilbert space Hi, and Jm is the couphng strength between the reference site 1 and 
all of the sites belonging to the k-th stratum with respect to 1. 

In order to use the spectral analysis methods, we write the hamiltonian (3-29) in terms of 
the adjacency matrix A of the network J {2m, m). To do so, first we note that for the Johnson 
network J(n, m) with adjacency matrix A one can show that 



l<i<j<n 



( \ 

m 



\ 2 / 



n — m 



+ 



v 



(3-30) 



where Pij is the permutation operator acting on sites i and j. In fact restriction of the operator 
J2i<i<j<n Pij on the m-particle space (space spanned by the states with m spin down) which 
has dimension C^, is written as the adjacency matrix A of the network as in the Eq. (3-30). 
For n — 2m, the Eq.(3-30) is written as 



E Pij + m{m - 1)1. 

l<i<j<2m 



(3-31) 



Then, by using the fact that 



(7j • aj = 2Pij - ^ ^ 2 ^ • = E ^'3-2 

l<i<j<2m l<i<j<2m. 



^ 2m^ 



\ 2 / 



I = A-^I, 
2 ' 



we obtain 



HG = J2JkPkiA). 

k=l 



(3-32) 



Now, for the purpose of the perfect transfer of a m- qubit GHZ state, we impose the con- 
straints that the amplitudes (0i|e'~*^*|0o) be zero for alH = 0, 1, m — 1 and (0m|e~'^*|0o) = 
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e*^, where 9 is an arbitrary phase. Therefore, these amphtudes must be evaluated. To do so, 
we use the stratification and spectral distribution associated with the networks J(2m, m) to 
write 

Let the spectral distribution of the graph is n{x) = Ylk=olkS{x — Xk) (see Eq. (2-18)). The 
Johnson network is a kind of network with a highly regular structure that has a nice algebraic 
description; For example, the eigenvalues of this network can be computed exactly (see for ex- 
ample the notes by Chris Godsil on association schemes [41] for the details of this calculation). 
Indeed, the eigenvalues of the adjacency matrix of the network J(2m, m) (that is x^s in ^i{x)) 
are given by 

Xk = — k{2m + 1 — A;), = 0, 1, . . . , m. 

Now, from the fact that for distance-regular graphs we have Ai = y/K7iPi{A) [27], {(J)., 
implies that 

m 

Y,-I^P,{xk)e-''^T=o-JiPi^-^) = 0, i = 0,l,...,m- 1 

fc=0 

Denoting e~**^i=o by r^jt, the above constraints are rewritten as follows 



iHt 



(3-33) 

|0o) = 



^Pi(xfc)%7fc = 0, i = 0, 1, ...,m - 1, 



fe=0 



k=0 



(3-34) 



From invertibility of the matrix Pj^ = Pi{xk) (see Ref. [2]) one can rewrite the Eq. (3-34) as 



Volo 



\ Vdld J 



' ^ 



(3-35) 



The above equation implies that 77^7^ for A; = 0, 1, m are the same as the entries in the last 
column of the matrix = WP* multiplied with the phase e'^, i.e., for the purpose of PST, 
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the following equations must be satisfied 

mik = 7fee-2**oSr=o'^'^'(-'=) = e'%WF')^^ , for A; = 0, 1, ..,m, (3-36) 

with W := diagi^Q, 7i, • • • , 7m)- 

One should notice that, the Eq.(3-36) can be rewritten as 

(Jo, Ji, . . . , Jd) = --^[^+(2Zo+/(0))7r, ^+(2Zi+/(l))7r, . . . , e+{2lD+f{D))T:]{W'P'), (3-37) 

or 

= + (2^. + /(J))^](W^P*).^, (3-38) 

where Ik for k — 0, 1, . . . , m are integers and f{k) is equal to or 1 (we have used the fact 
that 7fc and (M^P*)fcm ^^^^ for A; = 0, 1, . . . , m, and so we have 7^ = |(W^P*)feml)- '^^^ result 
(3-38) gives an exphcit formula for suitable couphng constants so that PST between the first 
node (|0o)) and the opposite one {\(t)m)) can be achieved. 

In the following we consider PST of the two qubit (the case m = 2) GHZ state l^^)^ = 
|(|00) -I- 1 11)) in details: Prom Eq. (2-23), for m = 2, the QD parameters are given by 

«! = 2, q;2 = 0; a;i = a;2 = 4, 

Then by using the recursion relations (2-13) and (2-16), we obtain 

Q^^\x) ^ - 2x - A, Qs{x)^x{x-A){x + 2), 

so that the stieltjes function is given by 

^ ( y_ QS\^) x^-2x-A 

Qs{x) x{x-i){x + 2y 

Then the corresponding spectral distribution is given by 

2 1 
^{x) = J2 liH^ - ^i) = r{3(^(2:) + 5{x - 4) + 2S{x + 2)}, 
1=0 " 
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which indicates that 



^ 70 ^ 



W 



71 

72 



^ 3 ^ 



1 

2 



In order to obtain the suitable couphng constants, we need also the eigenvalue matrix P with 
entries Pij — Pi{xj) — -^^==Qi{xj). By using the recursion relations (2-13), one can obtain 

Po{x) = 1, Pi{x) = f and Pa (a;) = \{x'^ - 2x - A) , so that 



2-1 
-111 



Now, from the result (3-38), we obtain 

69 + 6% 



TT 



71 
1% 



4 Conclusion 



By using spectral analysis methods and employing algebraic structures of Johnson networks 
J(2m, m) such as distance-regularity and stratification, a method for finding a suitable set 
of coupling constants in some particular spin Hamiltonians associated with the networks was 
given so that PST of m-qubit GHZ states between antipodes of the networks can be achieved. 
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